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Abstract 

Using the two dimensional XY — {S{0{y)) model as a test case, we show that analysis of the Fisher zeros of 
the canonical partition function can provide signatures of a transition in the Berezinskii-Kosterlitz-Thouless (BKT) 
universality class. Studying the internal border of zeros in the complex temperature plane, we found a scenario in 
complete agreement with theoretical expectations which allow one to uniquely classify a phase transition as in the 
BKT class of universality. We obtain Tbkt in excellent accordance with previous results. A careful analysis of the 
behavior of the zeros for both regions lHe(7’) < T bkt and lHe(7’) > T bkt in the thermodynamic limit show that 0fm(7’) 
goes to zero in the former case and is finite in the last one. 
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1. Introduction 

The Berezinskii-Kosterlitz-Thouless (BKT) transi¬ 
tion already has more than 40 years of history Ij]] and is 
still intriguing the scientific community. The nature of 
this transition is completely different from the common 
discontinuous (first order) or continuous (second order) 
phase transitions. Long range order does not exist and 
the two point correlation function has an algebraic de¬ 
cay at low temperature (T < T bkt) and an exponential 
decay for T > Tbkt 1211 ■ Here Tbkt is known as the 
BKT temperature, and a model displaying a BKT tran¬ 
sition has an entire line of critical points for T < Tbkt- 
In addition to that, the corresponding free energy, which 
is a C“ function, is not analytical in this region. Besides 
these striking features, the correlation function has a 
characteristic universal exponent decay ri(T bkt) = 1/4 
at the transition temperature. Its phenomenology relies 
on the belief that it is driven by a vortex-antivortex un¬ 
binding mechanism Hi . Another proposition that also 
describes the transition is based on a polymerization of 
domain walls J^]. Many systems, e.g., superfluid films. 
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Coulomb gases and crystal surface roughening, undergo 
transitions that can be classified as belonging to this uni¬ 
versality class 10]. More recently, Berker et al 10] found 
that a BKT transition can occur in a scale free network. 

Although the BKT transition is well known, the char¬ 
acterization of an unknown phase transition as being in 
the BKT universality class is not an easy task since there 
is no standard method to do so. The lack of a crite¬ 
rion capable of determining the nature of a transition 
beyond any reasonable doubt is a problem discussed by 
Bramwell and Holdsworth ]0]. They pointed out that 
to be able to see the transition the system under inves¬ 
tigation should be very large. They estimate that for “a 
system with atomic spacing of 3A the area should corre¬ 
spond to the size of a postage stamp”. From the analyt¬ 
ical point of view, the renormalization group approach 
is able to describe the main features associated with the 
transition. However, the approximations involved in the 
study of a given model may hinder the discovery of its 
real nature. 


In this paper we present an algorithm for the system¬ 
atic analysis of the Fisher zeros of the canonical par¬ 
tition function ]0, 10 HI] for the 2D XY model (with 
spin symmetry (9(3)) looking for possible signatures 
of the BKT transition. It was Fisher ifllll . in 1964, 
who proposed considering the partition function zeros 
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in the complex temperature plane to study phase tran¬ 
sitions. As known, the thermodynamic behavior of a 
given physical system is encoded by its partition func¬ 
tion, Z, and all thermodynamic quantities can be ob¬ 
tained as derivatives of the free energy, F - -ksT In Z. 
The basic assumption of the Fisher zeros approach to 
study phase transitions is that a system undergoes a 
phase transition at a given (real) temperature, Tc if 


T-1-1-1“ 


Z{Tc) - 0, reflecting the non-analyticity of F at Tc Ill2ll . 
Since in a BKT transition there is an entire line of crit¬ 
ical points for T < Tbkt, one should expect that a map 
of the Fisher zeros in the complex temperature plane ex¬ 
hibits an entire line of zeros in the real temperature axis 
for T <Tbkt, signaling the BKT behavior of the transi¬ 
tion, while for T > T bkt the zeros should not touch the 
real positive axis, emphasizing the analytical behavior 
of thermodynamic functions at high temperatures. Of 
course, these considerations apply only to the thermo¬ 
dynamic limit. We warn the reader that the Fisher zeros 
studied here should not be confused with the Yang-Lee 
zeros iSEl defined on the com plex fug acity plane in¬ 
stead of the temperature plane 113, 14, Q. In what 
follows we will analyze the Fisher zeros map for the 
XY-model. 


2. Model and Methods 

Here we study a “fruit fly” model of the BKT transi¬ 
tion |@], the classical two-dimensional XY-model on a 
square lattice, defined by 
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Figure 1: Fisher zeros map on the x = e“^'^ plane for the 2D classical 
XY-model in a 50 X 50 lattice. The inset shows a zoom on the inner 
region for five distinct simulations, represented by different symbols. 
The cusp, fj, is indicated. 


or the more reliable helicity modulus 11 Ei, quantities 
that do not assure the model is in the BKT universality 
class. 

In the Fisher zeros approach one starts with a discrete 
canonical partition function that can be written as Z = 
2 b giE) exp {-/3E), where g{E) stands for the density of 
states (DOS). For a continuous system one may perform 
a discretization of the DOS 12211 . This can be done by 
dividing the energy range into M bins of size s. We 
can organize the energies inside the interval [Eo,Em-i] 
by counting the energy of the bin as E„ = Eq + ne. 


By defining a variable x s e we get a discretized 
version. 


The sum runs over the nearest neighbors, J stands for 
the exchange coupling constant and S“ stands for the 
component a - (x, y, z) of the spin. The same Hamil¬ 
tonian also defines the Planar-Rotator (0(2)) model IB 
[B] , whose spins have only two components (one degree 
of freedom), and can also be viewed as an example be¬ 
longing to the BKT universality class. In spite of the 
lack of long-range order for the model H], there is a 
non-zero magnetization for any finite volume OlSll . re¬ 
sulting in a thermodynamic behavior very similar to that 
observed at continuous phase transitions. In fact, the be¬ 
havior can be easily confused with a second order phase 
transition or something else |B B 21]- Thus, distin¬ 
guishing between continuous and BKT behavior in fi¬ 
nite systems is difficult and may require system sizes 
beyond those feasible. Usually, the BKT temperature is 
estimated by using the Binder cumulant, the divergence 
of the magnetic susceptibility, the correlation functions 


M-l M-1 

Zo - e-^^» ^ g^x" = e-^^» Y]{x- Xj), (2) 

n=0 j=l 

where g„ = g{E„), and Xj is the /* zero of the polyno¬ 
mial, usually called Fisher zero. Once g„ is obtained, 
any reliable zero finder can be usedQ to find the x/s. As 
long as gn 6 Z is an analytical function and has no 
real positive zeros for finite systems. The zeros show up 
as complex conjugate pairs. The analysis of Fisher zeros 
in finite systems is done by considering the special set of 
zeros {x*(L) = a{L) -t- ib{L)] e Xj{L), where L is the lin¬ 
ear size of the system, called first or leading zeros. They 
have the following properties: b{L) —> 0 as L —> oo 


*In this work we decided to use: E.W. Weisstein, “Poly¬ 
nomial roots”, in MathWorld - A Wolfram Web Resource, 
http://mathworld.wolfram.com/PolynomialRoots.html 
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while lim /,^00 a{L) - a{oo), a constant value. The lead¬ 
ing zeros are very stable against statistical fluctuations, 
in contrast to non-leading zeros, and are, in gen eral, fea¬ 
tured in the map (see, for example. Ref. 12311 '). They 
are related to the transition temperature of the system in 
the thermodynamic limit as 1 /ksTc - - ln(a(oo')')/£ ll23n 
and their impact angle on the real positive axis is di¬ 
rectly related to the order of the phase transition ll24ll . 
Recent results also suggest that the pattern of zeros as a 
whole can be used to characterize the order of the tran¬ 
sition 


In a paper of 1983 Ytzykson et al. Il26l] . 
analysing the zeros distribution for small lattices, found 
that some properties of the Ising and gauge models seem 
to exhibit a universal behavior close to complex singu¬ 
larities. However, it is not clear how to extend their ar¬ 
guments to the BKT transition. 

In order to obtain the DOS of the XY-model we 
used the Replica Exchange Wang-Landau (REWL) 
method |27, 3 a parallel version of Wang- 

Landau (WE) sampling |3i, 32, 33], capable of sam¬ 
pling the entire configuration space efficiently in a sin¬ 
gle simulation. In this scheme the energy range is 
split into smaller overlapping windows. We considered 
eo = -1.97, cm-i - 0, and an overlap of 75%, where 
e - Ejl} stands for the energy per spin. Several differ¬ 
ent random walkers are allowed to run in each of these 
windows following the original WE schem^l We use 
regular square lattices with sizes ranging from L - \0 
up to 200. In this work we chose J - \,S - ks - 
and the lattice parameter a - 1. 


3. Results and discussion 

In Eig. [T]we show a typical zeros map of the imagi¬ 
nary and real parts of all x/s, for a 50 x 50 lattice. In 
a continuous phase transition a single leading zero is 
expected. However, we cannot identify a single point 
featured in comparison to others. The inset shows in 
different symbols the zeros for five different simulations 
for L - 50, in order to analyze statistical fluctuations. A 
leading zero cannot be identified in this picture. Instead, 
a cusp at is evident and the border line is quite stable 
against fluctuations. 

The size dependence of the internal border of the map 
pattern as a function of the lattice size is shown in Eig.|2] 
One should expect the internal border to coalesce into 
the real axis for 93e(x) < and L ^ 00 , in accordance 
with the existence of an entire line of critical points at 


^The considered flatness criteria is p = 0.7, the final In / value is 
10“®, and the acceptance ratio is 60%. 
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Figure 2: Zoom on the real positive semi-axis of the zeros maps in 
the X plane for L = 10 —» 100. The zeros on the internal border are 
highlighted. 


low temperatures. The line for fHe(x) > on the other 
hand should not touch the real axis. The cusp should, 
then, give Tbkt- 

In order to systematically investigate the finite size 
effects and confirm the above expectations, we opted to 
work in the complex temperature (T) plane instead of 


the complex x = e 


— p-s/ksT 


plane. This choice is justified 


by the fact that finite size scaling in terms of the tem¬ 
perature is known lll5n and that by doing so our results 
are almost independent on the chosen bin size, e. This 
last point is especially important when dealing with big¬ 
ger lattices. In this work we used e = 1 for L < 100 
and e = 3 for L - 200; we noted that small modifica¬ 
tions in the bin size did not change our results. To deter¬ 
mine the limits of the internal border we divided the real 
temperature axis into small bins of size centered at 
Tx - 9le(7’). The border line was identified by looking 
inside a given bin for the smallest value of Ty - 3m(7’) 
that appears. At least five different simulations were 
used for each lattice size. The resulting curves are pre¬ 
sented in Eig. |3] Eig. |4] shows Ty as a function of 
for some typical < Tbkt values. The solid lines are 
linear regressions showing that Ty ^ 0 for L —> 00 , i.e., 
the internal border coalesces with the real positive axis 
in the thermodynamic limit. Other scaling functions and 
exponents were also tested (not shown here), but they do 
not describe our data as well as this ansatz- 


To estimate the critical temperature we used the lo¬ 
cation of the cusp position, T^{L). Since the scaling 
function for the “pseudocritical” temperature is given 
by 115] r bkt{T) ~ [ln(L)]“^ we plotted T^{L) as a func¬ 
tion of [ln(L)]“^ in the Eig. A linear regression, 
discarding the point corresponding to L = 10, gives 
TBKT - 0.709(2), and discarding the points for L < 40 
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Figure 3: Internal border of zeros obtained by using the binning pro¬ 
cess for L = 10 ^ 200. Here ATx = 0.1 and the en'or bars represent 
statistical fluctuations. 



Figure 5: Typical Ty x L ^ behavior. Here Tx = 3.0 > Tbkt- Linear 
fit gives support to a power law behavior with exponent cj « 2.4. 



Figure 4: Finite size scaling analysis of the imaginary part of the in¬ 
ternal border according to the ansatz Ty ~ L~^. The lines represent a 
linear regression of the data, showing good agreement with the ansatz 
and the convergence to zero or small negative values for L ^ oo. We 
use IsTx = 0.1. 


gives Tbkt — 0-704(3), which agrees very well with 
previous results 1^, 18], 0.700(5) and 0.700(1) respec¬ 
tively. More precise results could be obtained by in¬ 
creasing the number of zeros in the map by reducing 
e, increasing the precision in the location of the cusp. 
However, the zeros finder task may become a problem 
for such a high degree polynomial (> 3 x 10"^). In any 
case, conventional methods to locate the BKT may be 
used together with this one to get even more precise re¬ 
sults for Tbkt- 


From Figs. |2] and [3 one can see that the curves di¬ 
verge from the positive real axis, in accordance with the 
expectation that for T > Tbkt the imaginary part of the 
zeros should remain finite. Moreover, since the free en¬ 
ergy has to be an analytic function at high temperatures, 
there can be no real positive zero of the partition func¬ 


tion at high temperatures. 



Figure 6: Finite size scaling of the cusp position, T,(L), according to 
the prediction for the BKT “pseudocritical” temperature Tbkt(L) ~ 
[ln(L)]“^. Discarding the point corresponding to L = 10 the BKT 
temperature is estimated by a linear regression as 0.709(2) (dashed 
blue line) and discarding L < 40, T bkt = 0.704(3) (solid red line). 


To analyse the behavior of Ty in the non-critical re¬ 
gion we have to be careful. First we have to consider 
the following. An analysis based in finite size scaling, 
as we did in the critical region, is meaningless since the 
basic assumption of any FS S is that the free energy be¬ 
haves as a homogeneous function. Following Ytzyk- 
son we identify three regions. (1) The region far from 
Tbkt where L » f ^ I, where finite size lattice ef¬ 
fects are not relevant. Here f is the correlation length 
which for the XY model diverges exponentially when 
Tbkt is approached from above, remaining infinite in 
the critical region. (2) The region T > Tbkt where 
L » f » 1. The finite lattice exhibits the same scal¬ 
ing behavior as the infinite system. And, (3) the region 
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where L x ^ » 1 for which are observed severe finite 
size effects. Regions (2) and (3) are scaling regions. 
With that in mind we content ourselves with an analysis 
in a region far enough from T bkt where we may expect 
that the border of the zeros map should converge fast to 
its asymptotic limit. For T > Tbkt the free energy is an 
analytical function of T so that it can be approximated 
by a polynomial in T with the coefficients depending 
on the lattice size in powers of 1/L. Farthest we are 
from T bkt faster will decay the higher coefficients of 
the polynomial, so that, it is reasonable to suppose that 
in this region Ty has a power law Ty(oa) + behav¬ 
ior. Asymptotic con verg ence is assured by the Brouwer 
fixed point theorem 13511 . Reasoning in this way we can 
plot Ty for the largest lattices of our simulations as a 
function of LT’^. By varying u) we can find the best value 
that adjusts a linear function to our data. This procedure 
is shown in Fig. |5] It is noteworthy the different “scal¬ 
ing” behavior of 3m{T) in both regions, T < Tbkt and 
T > T bkt- As should be expected the exponent at is not 
universal, depending on T. 

4. Conclusion 

In summary, we show that a method of investigat¬ 
ing the Fisher zeros of the partition function can iden¬ 
tify whether or not the model exhibits a Berezinskii- 
Kosterlitz-Thouless (BKT) phase transition. By study¬ 
ing the 2D XY-model we found a qualitative picture 
that is completely consistent with expectations for the 
BKT transition, i.e., the zeros map is consistent with 
the existence of an entire line of zeros in the real posi¬ 
tive axis in the thermodynamic limit, with different be¬ 
haviors for points below and above the transition tem¬ 
perature that could be used to signalize the BKT tran¬ 
sition. Moreover, the BKT transition temperature was 
successfully obtained by considering the location of the 
cusp that splits regions with different behaviors giving 
TBKT - 0.704(3), in excellent accordance with previous 
results llisl 13411 . Our quantitative analysis shows that 
for temperatures above T bkt the imaginary part of the 
zeros converge to finite values. This behavior is in ac¬ 
cordance with the fact that the free energy must be an 
analytical function over the real axis at high tempera¬ 
tures, which prevents the existence of any real positive 
zero at high T. 
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